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INTRODUCTION 


One  of  the  applications  of  studying  mean  stress  effects  in  fatigue  is 
residual  stress.  Often  in  the  manufacture  of  engineering  structures  residual 
stresses  are  induced  in  the  structure  either  intentionally  or  inadvertently. 
Examples  of  the  former  are  shot  p^ened  automotive  components,  swaged  holes  in 
aircraft  skins  or  autof rettaged  thick  cylinders.  Falling  into  the  latter 
category  are  welded  structures  or  certain  heat  treated  components  with  complex 
geometries.  When  such  a  component  is  subjected  to  fatigue,  its  life  can 
either  be  extended  or  decreased  depending  on  the  sense  of  the  residual  stress. 
Normally  compressive  stresses  result  in  enhanced  life  while  tensile  stresses 
usually  accelerate  failure.  Much  work  has  been  done  to  study  the  effects  of 
residual  stress  on  fatigue^  in  inert  environments.  If  the  additional  variable 
of  aggressive  environment  is  included,  things  become  even  more  complex.  To 
this  author's  knowledge,  no  organized  study  has  been  undertaken  to  address 
this  problem. 

Although  no  work  has  been  published  dealing  exclusively  with  the  residual 
stress  aspect  of  corrosion  fatigue  crack  growth,  much  work  has  been  conducted 
to  study  the  basic  phenomena.^  Fatigue  crack  growth  in  aggressive  environ¬ 
ments  is  influenced  by  many  factors  including  loading  frequency,  waveform, 
temperature,  and  mean  stress  or  stress  ratio  R  (R  =  ^min/^raax) •  these,  the 

R  effect  is  that  factor  which  is  related  to  residual  stresses.  This  aspect  of 

^Throop,  J.  F.,  Report  to  ASTM  Committee  E-9  on  the  Formation  of  Subcommittee 
E09.02  on  Residual  Stress  Effects  in  Fatigue,"  Philadelphia,  PA,  November 
1980. 

^Wei,  R.  P. ,  in  Fatigue  Mechanisms,  Jeffrey  T.  Fong,  ASTM  STP  675,  ASTM, 
Philadelphia,  PA7”1979,  p.  816. 
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environmentally  enhanced  fatigue  crack  growth  is  the  subject  of  this  report. 

The  aggressive  environment  system  used  to  study  these  effects  was  liquid 
mercury  and  a  low  alloy  high  strength  steel.  It  is  well  known  that  liquid 
mercury  very  strongly  attacks  steels  to  cause  liquid  metal  embrittlement 
(LME).3  Although  not  a  commonly  occurring  embrittlement  system  in  industry, 
the  system  is  very  similar  to  the  hydrogen-steel  system.  It  will  be  shown 
that  crack  growth  in  mercury  follows  the  same  shaped  curve  as  the 
hydrogen-steel  curve^  and  the  magnitude  of  crack  growth  rate  is  also  about  the 
same.  Using  mercury  thus  provides  a  good  analogy  to  a  practical  system 
without  the  expense  of  special  testing  chambers  and  equipment. 

We  will  show  the  results  of  tests  performed  at  various  values  of  constant 
R.  A  model  will  then  be  proposed  to  mathematically  fit  the  data.  Short¬ 
comings  of  the  model  will  be  discussed  suggesting  that  the  full  effect  of  the 
environmental  interaction  cannot  be  attained. 

EXPERIMENTAL  PROCEDURE 

The  basic  experiments  involved  in  the  testing  reported  herein  are  the 
measurement  of  crack  growth  rates  in  the  environment  under  controlled  R 
conditions.  It  was  decided  to  use  pure  ben4ing  samples,  since  these  can  be 
used  to  very  easily  produce  controlled  R  testing  conditions  and  are  easily 
manufactured. 


^Wei,  R.  P, ,  in  Fatigue  Mechanisms >  Jeffrey  T.  Fong,  ASTM  STP  675,  ASTM, 
Philadelphia,  PA,  1979,  p.  816. 

^Stoloff,  N.  S.,  in  Environment  Sensitive  Fracture  of  Engineering  Materials, 
Z.  A.  Foroulis,  Ed.,  The  Metallurgical  Society  of  AIME,  New  York,  1979,  p. 
486. 
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Since  we  were  concerned  with  only  the  R  value  effect,  the  initial 
experiments  were  performed  at  constant  waveform,  loading  frequency,  and 
temperature.  These  conditions  were  sine  wave  at  30  Hz  at  room  temperature. 
Five  R  conditions  were  tested,  R  =  ±  0.6,  ±0.3,  and  0.0.  In  a  cracked  body  R 
is  determined  as  the  ratio  of  stress  intensity  factors,  K,  rather  than  the 
ratio  of  stresses,  (R  =  K^in/Kniax)  •  The  subscripts  refer  to  the  minimum  and 
maximum  K  respectively  that  occurs  during  the  fatigue  cycle.  Figure  1  shows 
the  bend  sample  which  was  used  to  measure  the  fatigue  crack  growth  rates.  The 
K  solution  for  this  specimen  is  well  known  and  the  expression  for  K  as  a 
function  of  loading  conditions  and  crack  depth  is  given  as 5^ 


f(a/W) 


K  = 


4M 

B/w3/2 


f(a/W) 


3(a/W)^/^ 

2(l+2a/W)(l-a/W)372 


(1.99  -  a/W(l-a/W) 

X  (2.15  -  3.93  a/W  +  2.7(a/W)2)) 


(1) 


where  M  is  the  applied  moment,  and  a,  B,  and  W  are  the  crack  length 
thicknesses  shown  in  Figure  1. 


^Srawley,  J.  E. ,  International  J.  of  Fracture  Mechanics,  Vol.  12,  June  1976 
p.  475.  —  ^  . 
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Figure  1.  Schematic  Representation  of  the  Bend  Specimen  Used 
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All  of  the  variables  in  Eq.  (1)  were  fixed  once  the  testing  commenced 
except  for  the  crack  length,  a.  To  measure  crack  growth  in  aggressive 
®*^vlronments  several  methods  have  been  used.  In  this  study  we  used  the  crack 
mouth  opening  displacement  (CMOD)  method. 5  Using  the  technique,  the  CMOD  is 
continuously  measured.  With  all  other  variables  fixed,  as  the  crack  grows, 
the  CMOD  Increases.  Expressions  similar  in  form  to  Eq.  (1),  relating  CMOD  and 
the  other  loading  and  specimen  variables  for  many  specimens,  are  presently 
being  developed. ^  For  the  specimen  considered  here,  the  expression  is: 
a/W  =  6'(-0.980  +  5.1506'  -  4.2886'2  +  1.116' 3) 


6' 


15. 8M 

1  +  ( - ) 

(ACMOD)EBW 


1/2 


(2) 


where  E  is  Young's  modulus  of  the  material,  and  ACMOD  is  the  range  of  CMOD 
measured  in  the  experiments.  Equation  (2)  is  accurate  to  within  ±  one  percent 
for  any  value  of  6’ , 

The  specimens  were  prepared  such  that  LME  was  produced.  This  requires 
that  two  criteria  be  met:^  (1)  Intimate  contact  between  the  liquid  metal 
and  the  solid  metal,  and  (2)  tensile  stresses  at  the  crack  tip.  The  first 
criterion  Is  more  difficult  to  obtain  since  wetting  must  occur.  Steel  Is  very 


o 

Stoloff,  N.  S.,  in  Environment  Sensitive  Fracture  of  Engineering  Materials, 

Z.  A.  Foroulis,  Ed.,  The  Metallurgical  Society  of  AIME,  New  York.  1979,  p. 
486. 

^Yoder,  G.  R.  and  Crooker,  T.  W. ,  "The  Procedure  for  Precision  Measurement  and 
Analyzing  Fatigue  Crack  Growth  Using  Crack  Opening  for  Corrosion  Fatigue," 
presented  at  ASTM  E-9  Symposium,  Pittsburgh,  PA,  1979. 

Kapp,  J.  A.  and  Newman,  J.  C.  Jr.,  and  Gross,  B.,  in  progress. 
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difficult  to  wet  vd.th  liquid  mercury.  To  alleviate  this,  the  specimens  were 
first  plated  with  a  thin  copper  coating  which  can  be  wetted  very  easily  with 
mercury.  The  specimens  were  then  fatigue  loaded  to  break  the  copper  coating 
and  develop  fatigue  cracks  in  the  steel  beneath.  Since  freshly  created  crack 
surfaces  are  very  clean,  wetting  occurs  easily.  In  this  way  we  were  able  to 
produce  embrittlement  quite  readily.  The  steel  used  was  gun  steel  heat 
treated  to  a  yield  strength  of  about  1170  MPa. 


MATHEMATICAL  MODEL  OF  EMBRITTLED  CRACK  GROWTH 

Crack  growth  in  inert  environments  in  the  steel  tested  follows  the  well 
known  Paris  power  law:^ 


da 

--) 

dN 


=  CAK™ 


(3) 


inert 

where  da/dN)j[nert  the  crack  growth  rate,  AK  is  the  range  of  stress 
intensity  factor,  and  C  and  ra  are  empirical  constants.  It  has  been  suggested^ 
that  the  total  fatigue  crack  growth,  da/dN)cotal>  aggressive  environment 

can  be  modeled  as: 


da/dN)total  =  da/dN)inert  +  da/dN)envir  (4) 

where  da/dN)g^yj|^ is  the  contribution  of  the  aggressive  environment  which  can 
be  approximated  based  on  the  sustained  load  stress  corrosion  cracking  behavior 
of -the  material  in  the  environment  considered. 


^Paris,  P.  C.  and  Erdogan,  F. ,  Trans.  ASME,  Vol.  85,  1963,  p.  528. 
dWei,  R.  P.  and  Landes,  J.  D. ,  Mat.  Res.  Std. ,  Vol.  9,  1969,  p.  25. 
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Using  Eq,  (4)  we  are  able  to  develop  a  model  for  crack  growth  in  an 
aggressive  environment  which  incorporates  the  effects  of  R  values.  Consider 
the  idealized  sustained  loading  crack  growth  response  in  Figure  2.  If  the 
applied  K  is  below  some  critical  value,  Kxscc>  crack  growth  occurs.  When  K 
exceeds  KxscC>  crack  grows  at  a  constant  velocity,  v* ,  until  the  fracture 

toughness,  Kxc  is  reached  at  which  point  unstable  fracture  occurs.  If  we 
assume  that  <ia/dN)0^yXron  totally  accounted  for  by  this  property  of  the 
envi ronment ^me tal  system,  this  term  can  easily  be  determined  on  a  cycle  by 
cycle  basis.  For  one  cycle,  the  component  of  fatigue  crack  growth  is  given  as 
da 

”)envir  =  v*  x  (that  period  of  time  that  K  exceeds  Kxscc)  (5) 

To  determine  the  period  of  time  that  K  exceeds  KxscC  i^  s.ny  single  cycle, 
Figure  3  is  useful.  Referring  to  the  figure,  the  time  at  which  K  equals  KxscC 
is  t*.  The  amount  of  time  which  K  exceeds  KxscC  is  then  twice  the  amount 
(7T-t*)/2’TTf ,  where  f  is  the  frequency  in  Hz.  Thus  Eq.  (5)  becomes 

'  da  v*(Tr-t*) 

— )  - -  (6) 

dN  ,  TTf 

envir 


The  value  of  t*  must  be  calculated  to  determine  the  environmental 
contribution.  To  do  this  it  is  first  necessary  to  get  the  expression  for  K  as 
a  function  time  in  terms  of  the  variables  AK  and  R.  Referring  again  to  Figure 
3  we  find  that  for  one  cycle  at  one  Hz 
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Figure  2.  The  Idealized  Sustained  Load  Crack  Velocity  Property 
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Figure  3.  One  Sinusoidal  Fatigue  Cycle 


r 


AK 

K(t)  =»  ~  cos  2i\t  +  Kmean  (7) 

1 

^ean  * 

“  ^ax  "  ^in 

^in 
R  - - 

^ax 

It  can  be  easily  determined  that 

AK 

*  (I-R) 

RAK 

K.i„  -  --;j-  (9) 


AK  1+R 


Using  Eq.  (10),  Eq.  (7)  becomes 

AK  1+R 

K(t)  =  — ( - -  cos  2irt) 

2  1-R 


(10) 


(11) 


When  K(t)  =  KlSCC>  value  t* ,  corresponding  to  0  =  0*  =  2irt*,  making 

Eq.  (11) 

AK  1+R 

K(t*)  -  Kiscc  =  — ( - cos  0*)  (12) 

2  1-R 

Solving  for  0*,  Eq.  (12)  becomes 

( 1+R)  2Kiscc 

0*  =  arcos{ - - - }  (13) 


The  time  at  which  K  exceeds  Kxscc 


1  2(Tr-0*) 

time  =  -  X  - -  (14) 

f  2ti 
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Combining  Eqs.  (13)  and  (14) 


^  ,  1  (1+R)  2Kiscc 

time  that  K  exceeds  Kxscc  “  [l  “  “  arcos  { - —  - - ■}]  (15) 

TT  ( 1-R)  AK 


Combining  Eqs.  (15)  and  (6)  the  environmental  component  of  fatigue  crack 
growth  is 


da 

--) 

dN 

envir 


1  ( 1+R)  2Kiscc 

-arcos  [ - - - ]  ] 

TT  (1-R)  AK 


(16) 


There  are  some  conditions  which  are  applied  to  Eq.  (15)  to  make  it  valid. 
If  the  maximum  K  occurring  during  the  fatigue  cycle  is  less  than  KxscC>  there 
is  no  environmental  contribution.  Furthermore,  if  the  minimum  K  during  the 
fatigue  cycle  is  greater  than  KxscC>  there  is  always  an  environmental  effect. 
This  corresponds  to  t*  =  0.  Under  this  condition  the  environmental  component 
is  easily  obtained  from  Eq.  (6)  by  setting  t*  =  0.  Realizing  these  facts  and 
using  Eqs.  (8)  and  (9)  for  limits,  the  total  crack  growth  rate,  Eq.  (4)  can  be 
written  as: 


CAK™ 


AK 

(1-R) 


<  Kxscc 


da 

—  ^  J 

dN 


CAK®  +  --  [1 
f 


1  ( 1+R)  2Kxscc 

-arcos  [ - - - ]  ] 

TT  ( 1-R)  AK 


RAK 

(1-R) 


<  Kxscc 


AK 

< - 

(1-R) 


V* 

CAK®  +  — 
f 


KiSCC  < 


RAK 

(1-R) 


(17) 
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Using  this  equation  to  model  the  combined  effects  of  aggressive 
environment  and  stress  ratio  has  several  advantages.  If  it  indeed  applies, 
there  are  two  material  properties  which  need  be  found  to  explain  the  data,  the 
velocity,  v* ,  and  K^SCC*  This  being  the  case,  only  one  set  of  experiments,  in 
theory,  need  be  performed  in  order  to  develop  the  crack  growth  rate  property. 
Also,  using  an  equation  of  this  type,  crack  growth  predictions  can  be  made 
quite  easily  by  numerical  integration.  The  major  drawback  is  that  it  applies 
for  only  one  waveform,  sinusoidal. 

RESULTS  AND  DISCUSSION 

The  crack  growth  rate  results  for  five  different  constant  values  of  R  (0, 
±  0.3,  i  0.6)  are  shown  in  Figure  4.  The  straight  line  is  the  inert  environ¬ 
ment  property  of  the  material  tested,  which  has  been  previously  developed.^ 

The  curved  lines  are  fit  to  the  data,  using  Eq.  (17).  Several  comments  about 
the  data  can  be  made. 

It  is  apparent  from  the  actual  measurement  of  crack  growth  rate  that  the 
threshold  AK  for  accelerated  crack  growth  increases  with  decreasing  R  value. 
This  result  conforms  the  assumption  of  a  threshold  KxsCGj  which  is  basic  to 
the  formulation  of  the  crack  growth  model.  If  we  use  the  R  =*  0  as  a  base, 
Kiscc  would  be  that  value  of  AK  where  deviation  from  the  inert  environment 
behavior  is  observed.  For  the  data  in  Figure  4,  this  value  is  about  33  MPa/ra. 
Eq.  (17)  predicts  that  accelerated  crack  growth  would  occur  for  non-zero 

^Underwood,  J.  H.  and  Throop,  J.  F.,  Part-Through  Crack  Fatigue  Life 
Prediction,  ASTM  STP  687,  J.  B.  Chang,  Ed .  ,  "^TM,  "Philadelph^^  1979,  p. 

195. 
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Figure  4.  Comparison  of  the  Measured  Results  with  the  Uncorrected  Model 
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R  values  when  AK  =*  33(1-R).  For  R  of  +0.6,  +0.3,  and  -0.6,  the  critical  AK 
values  are  calculated  as  13,  23,  33,  43,  and  53  MPaV^,  respectively.  The 
actual  critical  AK  values  for  the  respective  R  values  tested  were  15,  25,  33, 
44,  and  53  MPa/m.  This  is  excellent  agreement. 

The  crack  growth  behavior  beyond  the  critical  AK  value,  phenomenologi¬ 
cally,  can  be  explained  as  follows.  As  AK  increases,  da/dN  Increases  very 
rapidly  over  a  short  range  of  AK.  Increasing  AK  further,  the  Increase  in 
da/dN  is  much  more  gradual,  until  at  some  very  high  AK  value,  in  this  case 
about  100  MPa/m,  there  seems  to  be  little  effect  of  the  aggressive  environ¬ 
ment. 

The  general  shape  of  the  observed  behavior  is  well  characterized  by  the 
proposed  model.  Eq.  (17)  shows  that  there  should  be  a  limit  to  environmental 
effect,  when  AK  >  KxscC^ while  the  inert,  Paris  power  law  component 
always  Increases  with  AK.  This  indicates  that  there  could  be  some  value  of 
AK,  at  which  the  environmental  component  of  embrittlement  becomes 
insignificant.  The  model  predicts  this  limiting  behavior  quite  well  for 
negative  values  of  R  tested. 

When  R  is  positive,  the  maximum  crack  growth  rate  at  high  AK  values  is 
overestimated,  using  Eq.  (17).  In  this  case,  it  is  possible  for  K  to  exceed 
Kiscc  entire  fatigue  cycle.  The  model  proposed  by  Eq.  (17)  suggests 

under  these  conditions  that  the  crack  is  always  growing  at  the  velocity  v* 

from  the  environment,  in  addition  to  the  crack  growth  from  the  inert  fatigue 

conditions.  The  actual  experimental  results  shown  in  Figure  4  suggest  that 

there  is  a  limit  to  the  amount  of  environmental  interaction  which  occurs 
regardless  of  R  value.  Please  note  how  all  the  data  converge  to  a  crack 
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growth  rate  of  about  10~^  m/cycle  at  AK  of  about  100  MPaVm.  Equation  (17) 
predicts  that  at  AK  of  100  MPaV^  the  crack  growth  rate  approaches  1.6  x  10~^ 
m/cycle  for  either  positive  R  value  testing  condition.  This  suggests  that 
only  about  63  percent  of  the  full  environmental  effect  was  actually  attained 
in  the  experiment. 


The  overestimation  of  the  actual  behavior  also  occurs  when  a  partial 


environmental  effect  is  expected.  For  example,  consider  the  R  =  +0.3 
condition  at  AK  of  about  50  MPa/m.  The  total  crack  growth  rate,  given  by  Eq. 
(16)  is  about  7xl0~^  m/cycle.  The  inert  component  is  8.15xl0~^  m/cycle  and 
the  environmental  component  is  6.185x10“^  m/cycle.  If  the  environmental 
component  is  reduced  to  63  percent  of  the  predicted  value,  the  total  crack 
growth  rate  would  be  4.7x10  ^  m/cycle.  In  the  experiment,  the  crack  growth 
rate  measured  was  about  4.4x10“^  m/cycle,  which  agrees  quite  well  with  the 
predicted  value.  Therefore,  a  better  fit  to  the  actual  data  is  suggested  as 
follows : 

r 

CAK™  AK  <  (l-R)KiscC 


da  Av*  1  ( 1+R)  2Kiscc 

/  CAK®  + - [1  -  -arcos{ - }] 

dN  1  f  IT  ( 1-R)  AK 

,  ,  ( 1-R) 

(l-R)KxscC  AK  < - KiSCC 

R 


(18) 


Av* 

CAK®  +  - 

f 


(1-R) 

AK  > - RiSCC 
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where 


r 


1.0  R  <  0 


A 


0.63  R  <  0 
C  =  6.52  X  10-12 

m  =»  3 

f  =  30  Hz 

l^ISCC  *  MPa/m 

Using  Eq.  (18),  the  actual  observed  data  is  compared  with  theory  in  Figure  5 
showing  much  better  agreement. 

The  reasons  for  this  type  of  behavior  at  positive  R  values  are  not  immedi¬ 
ately  apparent.  A  possible  explanation  may  be  an  inability  of  the  liquid 
metal  to  access  the  crack  tip  for  embrittlement  to  take  place.  Gordon^^  has 
suggested  that  the  rate  of  crack  growth  in  liquid  metal  embrittlement  may  be 
limited  by  the  flow  of  liquid  metal  in  capillary  created  by  the  newly  created 
fracture  surfaces.  Using  fluid  mechanics,  an  equation  was  derived  which 
related  the  fluid  velocity  (Vf)  to  the  separation  of  the  fracture  surfaces, 
approximated  as  the  crack  tip  opening  displacement  (CTOD) : 

Vf  =  (CTOD)^  •  f (solid  metal-liquid  properties,  crack  length)  (19) 
Under  static  loading  condition,  CTOD  is  a  function  of  the  yield  strength 
of  the  solid  metal  and  applied  K.  If,  under  the  cyclic  fatigue  conditions 
studied  here  we  consider  the  CTOD  to  be  an  RMS  value  for  a  sine  wave,  the 
CTODrms  is  about  0.707  of  the  maximum  CTOD.  Therefore,  the  fluid  flow 
according  to  Eq.  (19)  under  cyclic  loading  would  be  0.5  of  the  CTOD  under 


lOCordon,  P. ,  Met.  Trans.  A.,  Vol.  9A,  1970,  p.  614. 
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do/dN  (m/cycle) 


AK  (MPaym) 


Figure  5.  Comparison  of  the  Measured  Results  with  the  Modified  Model 
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static  conditions.  If  it  is  assumed  that  the  fluid  velocity  limits  the  crack 
growth  rate,  the  above  analysis  suggests  that  only  about  50  percent  of  the 
full  environmental  effect  could  be  attained.  From  the  experimental  data,  it 
was  found  that  about  63  percent  of  the  full  environmental  effect  was  observed. 
This  comparison  suggests  that  the  fluctuating  nature  of  the  crack  tip  opening 
displacement  may  account  for  the  need  for  correction  factors  in  the  model  when 
R  is  positive.  It  also  may  suggest  that  this  effect  may  be  needed  when  R  is 
negative  because  the  waveform  is  also  sinusoidal  in  these  cases  and  the  CTOD 
is  certainly  changing.  In  such  instances  only  a  portion  of  the  fatigue  cycle 
is  above  the  Kjgcc,  and  the  environmental  interaction  can  only  occur  during 
that  portion  of  the  cycle  when  K  is  greater  than  Kiscq.  Under  these 
circumstances,  the  CTOD  does  not  change  as  much  as  when  the  entire  cycle  is 

above  Kjgcc.  Thus,  the  correction  from  the  model  should  be  less  in  the  cases 
when  R  is  negative. 

CONCLUSIONS 

The  effect  of  mean  stress  or  R  ratio  on  fatigue  crack  growth  in  an 
aggressive  environment  has  been  measured.  The  crack  growth  rates  can  be 
accounted  for  very  effectively  by  a  simple  superposition  model  when  R  is  zero 
or  negative.  A  correction  must  be  made  to  the  model  when  R  is  positive.  The 
justification  for  the  correction  factor  is  the  restriction  of  the 
environmental  access  to  the  crack  tip  from  cyclic  loading. 
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